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0^ | We have studied the penetration depth for a model of Y Ba2Cu^Oi involving pairing both in 

the CuOi planes and in the CuO chains. In this model pairing in the planes is due to an attractive 
' interaction, while Coulomb repulsion induces in the chains an order parameter with opposite sign. 

1) , Due to the anticrossing produced by hybridization between planes and chains, one obtains a d-wave 

■ like order parameter which changes sign on a single sheet of the Fermi surface and has nodes in 
the gap. We find that our model accounts quite well for the anisotropy of the penetration depth 
and for the absolute values. We reproduce fairly well the whole temperature dependence for both 
the a and the b directions, including the linear dependence at low temperature. We use a set 
of parameters which are all quite reasonable physically. Our results for the c direction are also 

■ satisfactory, although the situation is less clear both experimentally and theoretically. 
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The debate about the mechanism of high T c superconductivity has seen recently marked progress by focusing on the 
symmetry of the superconducting order parameter. This question is intimately related to the mechanism but lends 
itself much more easily to experimental answers. A number of recent experiments, mostly on Y BazCu-sO-i (YBCO), 
have shown that the order parameter displays lines of nodes on the Fermi surface. The most spectacular ones are 
the phase sensitive experiments [0J^] which provide striking evidence for a change of sign of the order parameter 
at the Fermi surface. However these experiments are sensitive to the quality of the interfaces which are not under 
perfect control. It is therefore important that other experiments, including tunnelling Raman scattering p| and 
penetration depth H, have shown the existence of low energy excited states. Although these experiments are not 
phase sensitive and could be explained by a strongly anisotropic s-wave order parameter, this interpretation seems 
fairly unlikely in view of the linear slope of the penetration depth at low temperature which is particularly striking 
and has been observed both in cristals || and in films 



A natural interpretation of the existence of nodes for the order parameter is the spin fluctuation mechanism J7[ which 
predicted a d x 2_ y 2 symmetry. But this hypothesis remains controversial because it meets with experimental as well as 
theoretical difficulties. The most prominent experimental stumbling block is the observation of a sizeable Josephson 
current Q in c-axis tunnelling between Y Ba^Cu-^O-j and Pb which contradicts a pure d-wave symmetry, and which 
is hard to interpret as being caused by the small orthorhombic distorsion because of the large size of the effect. Other 
experiments display also a strong ab plane anisotropy for the superconducting properties, such as the magnetic field 
dependence of the specific heat anomaly at T c of the parent compound LuBa 2 Cu^,Oy . Most strikingly the penetration 
depth has a marked anisotropy H since, at zero temperature, the ratio between the superfluid densities A^ 2 (0) / A~ 2 (0) 
is 2.3, in good agreement with the anisotropy of the normal state resistivity and of the square of the plasma frequency 
In Y Ba-iCu^O^, which has an orthorhombic distorsion of 0.8% (that is somewhat less than Y Ba-iCu^O-j which 
has 1.6%) and where there are twice as many chains as in YBaiCu^O-j, the superfluid density anisotropy Q is as high 
as 6.2 . In view of this evolution with increasing the number of chains, it is hard to escape the conclusion that the 
chains are conducting and that it is their direct contribution which is responsible for the superfluid density anisotropy. 
Now the linear term in X^ 2 (T) at low T is typically 3 times the linear term in A~ 2 (T). If we attribute the linear term 
to the existence of nodes in the gap, this implies, roughly speaking, that there must be something like nodes on the 
chains. This does not fit easily in the spin fluctuation picture and is rather a strong indication that chains must be 
explicitely taken into account for a proper description of the superconducting properties, even if it is believed that 
they are unlikely to be the place for the dominant driving force toward superconductivity. 

Another conspicuous problem of the spin fluctuation interpretation is the weak sensitivity of the critical temperature 
of Y Ba^Cu^O'j to the presence of impurities. Indeed standard impurities produce in d-wave superconductors p0[ an 
effect analogous to pair-breaking by magnetic impurities in standard s-wave superconductors with a rapid decrease 
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of the critical temperature with increasing impurity concentration following the Abrikosov - Gorkov law. However, 
except for Zn which is likely to have a magnetic environment once in an YBCO matrix, impurities seem to have a 
rather weak effect on the T c of Y Ba^Cu-sOi . Ion or electron irradiation have also shown a similarly weak sensitivity 
of the critical temperature. 

In order to account for all these apparently contradictory experiments we have proposed recently a model |Tl| where 
an order parameter with nodes is produced, not by in-plane repulsive interaction as in the spin fluctuation mechanism, 
but rather by a repulsive plane-chain pairing interaction which we attribute to mostly unscreened Coulomb repulsion. 
As a result the plane band and the chain band have order parameters with opposite sign. Since it is well known that 
these plane and chain bands hybridize because electrons are physically allowed to jump from planes to chains, this 
change of sign gives automatically rise to nodes of the gap in the region of the Brillouin zone where they anticross. 
In our model, although we do not exclude some pairing interaction coming from the chains, superconductivity arises 
primarily from an attractive interaction within the planes, due for example to phonon exchange. We can say that 
our order parameter is d-wave like because it is qualitatively similar |ll|] to the d x 2_ y 2 model : on each sheet of the 
Fermi surface there are four nodes. On the other hand, since the chains play an essential role in our model, it is 
meaningless to consider an approximate tetragonal symmetry. Rather we have to consider only the orthorhombic 
symmetry, under which our order parameter is completely invariant (there is no symmetry breaking). Note that the 
d x 2_ y 2 order parameter is also invariant under this symmetry. 

Since the order parameter of our model has nodes, it accounts naturally for all the experiments providing evidence 
for its change of sign and for the existence of low energy excitations. On the other hand, since our model is far from 
having the d x 2_ y 2 symmetry because of the important role of the chains, it provides a simple explanation for the 
existence of the Josephson effect in c-axis YBCO — Pb junctions and also naturally for the important anisotropy in 
the ab plane mentionned above. The existence of two weakly coupled bands (plane and chain) in our model leads to 
the possibility of a weak sensitivity of the critical temperature to impurities, in agreement with experiment, as we 
have shown very recently 113] . We note that our model has similarities with the two band S-N model, introduced 
by Abrikosov and Klemm p3[ to account for Raman scattering data. It has also common features with single band 
models, where competing attractive and repulsive ( Coulomb ) interactions lead to a change of sign of the order 
parameter within the band, such as the ones proposed by Abrikosov H] and by Santi et al. p5| . 

In this paper we will explore in detail the consequence of our model for the penetration depth. We note that, 
since it is a thermodynamic quantity, the penetration depth is less likely to be perturbed by extrinsic defects than a 
dynamical quantity. Moreover since one explores the sample over typically 1000 A, it is a bulk quantity. This is in 
contrast with surface experiments, such as photoemission or tunnelling, where only a few atomic layers are involved 
and which could give skewed informations if the surface happens to be somewhat different from the bulk for some 
reason. All these features make the study of the penetration depth particularly interesting. Since our model gives 
naturally nodes in the gap, it is clear that it gives rise qualitatively to a linear dependence of the penetration depth 
A(T) at low temperature. However it is not obvious that our model can account quantitatively for the size of the 
slope nor for the experimentally observed temperature dependence and its anisotropy. The purpose of the present 
paper is to explore these points in details. Here we will consider only Y Ba2Cu^O'j , since this compound is the best 
controlled and has been the most explored experimentally. Possible extension of our model to other compounds is left 
for future work. 

Experiments on the penetration depth in YBCO, whether on cristals or on films 0, are now in reasonable 

agreement. In particular they find a linear behaviour at low temperature with a typical slope of 4.3 A/K for an average 
penetration depth A a fc. More recently the anisotropy of the penetration depth has also been measured in untwinned 
single cristals |l7) . Since these last experiments provide stronger constraints on theory, we will mostly take them for 
comparison with experimental results ( note that the compound is more precisely Y BCI2CU3O6.95 but we will consider 
that it does not differ appreciably from Y BaiCu^O-j ). The T = values of the penetration depths are more difficult 
to obtain and we will make use of the results of Basov et al. ]9| for comparison. However experimental agreement is 
not complete and there might still be some progress in the future on this side, with perhaps some surprises Jl8| . 

II. THE MODEL 

We first introduce the essential ingredients of our model. The non interacting part of the Hamiltonian for the 
electrons in planes and chains is given pl| by : 
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H = ^ £ kctck + Y £ 'k d k d k + Y t k c k d k + h.C. (I) 
k k k 

where c k and d k create electrons in the plane and in the chain band respectively. The first term corresponds to 
an isolated plane while the second one describes isolated chains with respective dispersion relations tk and e' k . The 
last term describes hopping between planes and chains. Actually since YBCO is made of stacks containing two 
CuOi planes and one CuO-chains plane, this independent electron Hamiltonian does not correspond precisely to the 
situation found in YBCO. A more realistic description is obtained with the Hamiltonian: 

H = E £C i.n C i-n + E £C Jn C 2,« + E £p(c+„C 2 ,„ + h.C.) 
n n n 

+ E £ ' d t d n + E *c(c^„d n + h.c.) + J2 *c(c^ n d n+ i + h.c.) 

n n n 

Here all quantities are understood to depend on k x and k y and summations run also over k x . y . Indices I and 2 number 
the CuOi planes and the index n numbers the stacks. Introducing the even and odd plane band operators c± by c\.i 
= ( c + ± c_ )/ y/2 , and taking the Fourier transform in the z direction, we have: 



H o = E( £ + h) c % k c +,k + E( e ~ hY- k c ~,k+ 
k k 

+ E £ '4 d k + E M< fe 4 + h.c.) + J2t_ (c± k d k + h.c.)(3) 



(3) 



k 



Now a dependence and summation on k z is understood, and we have set t + = t c ^2 cos ( k z c / 2 ) and t _ = t c J2 
sin ( k z c / 2 ) , with c being the size of the unit cell along the z direction. From band structure calculations |19|-|2l|l 
the chain band Fermi surface crosses the odd plane band Fermi surface, whereas there is no crossing with the even 
plane band. Therefore the even plane band does not play an interesting physical role in our model. If it is omitted 
this leads back to the Hamiltonian Eq.(l) where the plane band is actually the odd plane band and t* = t c \/2 sin 
( k z c / 2 ) . In order to simplify the discussion we will forget the even plane band. However when we will come to 
realistic calculations we will take into account its contribution. 

With respect to pairing interactions, we take an attractive pairing in the plane while we assume a repulsive pairing 
between plane and chain. On the other hand we have no precise indication on the pairing interaction within the 
chains ( a simple and reasonable hypothesis would be to take the coupling constant equal to zero, but this is clearly 
somewhat arbitrary ). This leads to the effective interaction Hamiltonian ( with g and K positive ) : 

Hmt = ~ff y^Cfc/clfc/C-fcCfe + K^d^dt^c-kCk + h.c. - g f y]d£,dt k ,d-kdk (4) 

and the following mean field Hamiltonian : 

H = F + A^c+ci fc + A>J2dtd± k +h.c. (5) 

k k 

where the order parameter (A, A') satisfies : 

A=-g^<c_ k c k >+K^<d_ k d k > (6) 

k k 



A' = K < c_ k c k > -g' Y < d -kd k > (7) 

k k 

Because of hybridization, the plane and chain operators c k and dk do not correspond to the eigenstates of Hq and one 
has to perform a unitary transformation to diagonalize Hq . The energies e± (k) are given by : 

2e±(k) - e(k) +e'(k) ± [(e(k) - e'(k)) 2 + 4t 2 (k)] 1/2 (8) 

When this transformation is carried out in the pairing part of the Hamiltonian one obtains interband pairing terms 
coupling the hybridized bands. However we will consider that t & is large enough so that these bands are well separated. 
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Specifically this means that we assume A , A' << t k ( the parameters we use below do indeed satisfy this constraint, 
except for a small range around k z = which is quantitatively unimportant ). In this case one can easily see 
that these interband pairing terms are negligible. Therefore the transformation gives rise to a fully band-diagonalized 
Hamiltonian, with the order parameter in each band given explicitely by A±(k) = ( A e'(k) + A' e(k) )/( e(k) + 
e'(k)) ( k being at the Fermi surface of the + or - band ). The essential feature of our model is that, starting from 
isotropic interactions in planes and chains, we obtain a specific anisotropy for the order parameter ( A + (k) , A_(k) 
). In particular we obtain an order parameter which changes sign and has nodes at the Fermi surface since we have 
managed to have A and A' with opposite signs, while e(k) /(e(k) + e'(k)) goes essentially from to 1 when we move 
at the Fermi surface of a given band. Therefore we have an order parameter which is d-wave like, although we have 
assumed an attractive pairing in the planes. 

In order to obtain the temperature dependence of the penetration depth, we need the temperature dependence of 
the order parameter. This should be obtained as the solution of our two band model. As it is well known from the 
early work of Suhl ct al. j2^], the resulting temperature dependence can be appreciably different from the standard 
BCS result. However we will not require in our calculations a ratio |A'| / A between the chain and plane order 
parameters very different from unity ( it will be typically 0.5 ). We have checked that, in the range of coupling 
constants which is of interest for us, within a weak coupling calculation, the departure from the BCS temperature 
dependence is rather small. 

Another complication is that, to obtain the order parameter and the penetration depth, we should perform a strong 
coupling calculation by solving Eliashbcrg equations. This is necessary to be consistent with the high value of the 
critical temperature ( and indeed recent calculations of T c within the spin fluctuation model are also performed in 
the strong coupling regime p3fl ). This is also consistent with the rather high values of 2 A / T c whenever it has 
been measured experimentally, by tunnelling or Raman scattering for example. And indeed we will make use below 
of values for 2 A / T c somewhat larger than the BCS one. Such calculations do not make in principle any problem 
and they would be appropriate if we wanted precise quantitative results in order to fit experiment. However they 
have the disadvantage of introducing additional parameters - these would be the parameters necessary to describe 
the phonon spectrum |24|) , in addition to the coupling constants ( including intraband Coulomb repulsion ). However 
our purpose in the present paper is not to perform refined calculations, but rather to show that our model can give 
a reasonable quantitative account for experiments without looking for perfect agreement ( although as we will see we 
come already rather close to it ). In addition we will find other sources of uncertainty. In this spirit we will take 
the following simplifying approach, as it is often done in the literature. We will perform a weak coupling calculation 
for the penetration depth and assume the weak coupling BCS result for the temperature dependence of the order 
parameter. We will take strong coupling effects into account only by a renormalization of the ratio 2 A / T c . It is 
known that, for the moderate coupling constants we are interested in, strong coupling calculations do not give much 
different results. We will also invoke strong coupling effects when we will look for absolute values of the penetration 
depth since, because of mass renormalization, they will modify the bare results. The advantage of this approach is 
that it leaves us with only two adjustable parameters, namely the ratio A' / A between the chain and plane order 
parameters ( we assume each one isotropic for simplicity as in our original model ) and the ratio 2 A / T c . Moreover 
for both of these parameters we do not have physically much freedom. The price is a quantitative inaccuracy which 
we believe is rather small. 



III. IN-PLANE PENETRATION DEPTHS, LOW TEMPERATURE BEHAVIOUR 

Our calculations are made for clean superconductors. This limit seems appropriate for the experiments we consider, 
and more generally it seems to be the relevant one for high Tc superconductors because of their short coherence length. 
It would naturally be of interest to study the effect of impurities on the low temperature behaviour of the penetration 
depth, since this can help to discriminate between various theories. However since the physics is in this case controlled 
by the vicinity of the nodes of the gap, we do not expect much difference, at least qualiltatively, with what is found 
for the standard d-wave model |p5| . On the other hand it is not clear that impurity effects can be entirely neglected 
in experiments, in particular for the penetration depth along the c-axis (we will take up this point again when we will 
look at this penetration depth ). Nevertheless we do not consider them in the present paper. 

As we have already mentionned, in order to account for the experimental results we will need a fairly sizeable value 
of the hybridization parameter t% , which will be much larger than a typical order parameter. In this case the two 
hybridized bands are well separated and we have to consider a multiband situation. In such a case, since we can not 
use a free electron description, one may worry about the contribution of interband terms to the paramagnetic current 
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as in Ref. [^6). However one can check that the superfluid density is zero at T c as it should ( this results directly 
from the "effective mass theorem" |27j ), that is the paramagnetic current is exactly cancelled by the diamagnetic 
one. Below T c the modification of the interband terms due to the superconducting condensation is small provided 
the separation between the bands ( that is the difference between the energies of one-electron eigenstates with same 
wavevector ) is large compared to the gap. More precisely, if D is a typical band separation and A is a typical gap, 
the correction is of order ( A / D ) 2 , and should be neglected for consistency in our case ( D is of order of the Fermi 
energy Ef except in the anticrossing region where it is or order tk , but since this region is small, the correction is of 
order (A / Ep) 2 anyway ). 

The expression of the penetration depth in this case is just the sum of the contributions from each band. This leads 
to the standard result along the i = x or y axis: 

K\T) [ d*k ^ Al(T) 

-^-5 = / ^o(e±)v i ± 2nT ) -p- (9) 

where the summation is over both the + and the - bands, and is limited to the Fermi surfaces by the Dirac function 
5(e±). The gaps A± , the energies e± and the velocities v^± are naturally all k - dependent with the velocities 
u,- l ±(k)= de±(k)/dki . The discrete summation in Eq.(9) is over the Matsubara frequencies w n = (2n+l)7rT. 

We will first consider the low temperature behaviour of the penetration depth. Indeed on one hand the most striking 
feature of the recent experimental results is the linearity of the low T dependence. On the other hand we can obtain 
in this regime explicit expressions in our model with a limited number of assumptions, because this dependence 
is controlled by the nodes of the gap which are located in the anticrossing region. We will obtain an analytical 
expression for this low temperature behaviour in the limit of small hybridization. With this aim we note that there is 
no linear term at low temperature if hybridization is zero. Hence if we subtract from Eq.(9) its counterpart with the 
hybridization parameter tk set to zero, the low T behaviour is unchanged. Then we are left with summations over 
quantities which are essentially nonzero only in the anticrossing region. If we assume tk small enough, this region 
is small and we can perform a first order tk expansion ( since the unhybridized term does not contribute to the low 
T, we can forget it in the rest of the calculation ). In this anticrossing region the Fermi velocities <9e(k)/9k and 
<9e'(k)/9k of the plane and chain bands can be taken as constants, equal to their value v p and v c at the crossing of 
the unhybridized Fermi surfaces. 

It is more convenient to make a change of variables ( at fixed k z ) and take x = ( e - e')/tk and y = ( e + e')/2tk as 
new variables instead of k x and k y . The Jacobian of the transformation is J = d(x, y) / d(k x , k y ) = | v p x v c | / i 2 . . 
The y integration is then easily performed. One can then see, by changing x into -x, that the contribution of the two 
bands + and - are identical, and we are thus left with the calculation on the - band. The low T behaviour comes from 
the vicinity of the node x of A_ = (1/2) [ ( A + A' ) - ( A - A')x/(x 2 +4 f/ 2 ). This is given by x = ( A - | A'| ) 
/ ( A I A'| ) 1/2 . Around this node, A_ = -2{x - x ) ( A |A'| f' 2 / ( A + | A'| ) 2 . In performing the x integration, 
we obtain the dominant T dependence by evaluating Vi t - at the node xq . Since the gap A_ is now the only varying 
quantity, we are now led to a situation which is essentially identical to a d-wave calculation : after taking A_/T as 
new variable, we are left with a purely numerical integral which turns out to be equal to tt ln2 . Taking into account 
the two bands ±, we obtain for the low temperature dependence : 

Ar>) -^ ( T) = i^lT<™! (&) (10) 
* n 2 c (A|A'|) 3/2 V Jt k/ 

where c = 11.7 A is the size of the cell along the c axis, t is the average of |ifc | over k z and the value of (vf _ / J tk 2 
)o at the node is given by : 



2 \ /A , I A II \2 



V l.- 



Jt 



(Av c ,i + |A'|vp,i 
lj n (A+|A'|) 2 |v p xv c 



(11) 



This dimensionless quantity is typically of order unity. Expressing numerically all known physical constants, this can 
be rewritten as : 



v 3, n ,TT c (1 + |A'|/A) 



2 / „2 



A i( T) - A, ! ) 1.7(^(11) — — ", ; "| %r 1 (I 2 ) 



Tc A f| A '|/A) 3/2 V Jt k 



5 



where all the penetration lengths are now expressed in units of 1000 A and the hybridization parameter t is expressed 
in eV. We want now to compare this result with experiments |H| on YBC12CU3O7 which gives a slope of 4.7 A/K in 
the along the a axis, and 3.6 A/K along the b axis. 



A first point to make is that our result contains the zero temperature value Ai(0) of the penetration depth. We 
can naturally perform a calculation for this quantity, as we will do below. However the necessary ingredients to 
obtain numerical results, namely band structure informations and mostly rcnormalization due to the interactions, are 
not precisely known. Moreover Ai(0) is found by performing an averaging over the whole Fermi surface. Hence its 
physical content is essentially unrelated to the other ingredients of Eq.(12) for which only the vicinity of the nodes 
is relevant : they give a very local information while Ai(0) contains a global one. Therefore we might just take Ai(0) 
from experiment itself, which gives H the currently accepted values of A a (0) = 1600 A and A(,(0) = 1000 A. This is 
all the more reasonable as one sees that our result is fairly sensitive to the value of Ai(0) since it comes in the third 
power. 

Even taking A^(0) from experiment we have many parameters in our result which are not so well controlled. We 
can first consider the anisotropy of the linear term, since most of the factors in Eq.(12) are common to the two axes. 
We have : 



Since the chains are along the y direction, it is quite reasonable to assume that the component v cx of the chain 
electron velocity along x is essentially zero . Setting v c ^ y — v c , we obtain the following numerical condition between 
our parameters in order to agree with experiment : A / | A' | + v p . y / v c — 1.77 v ViX / v c . This relation is easily 
satisfied with gap ratios and velocity ratios of order unity, which is just what one would expect. It is therefore not a 
strong constraint. Just to give an example, the simple minded tight binding dispersion relations ek = — 2 t\ ( cos k x 
+ cos ky ) — jU and e' k = — 2 t\ cos k y — // for the plane and the chain bands give A / | A' | = (1.77 sin k x / sin k y 
— 1 ) t\ 1 1[ . This is satisfied with t\ w^a crossing located |^0) at k x ~ n/2 and k y w ir/4 , and | A' | / A « 0.7 . 

We consider now if we can account for the absolute value of the slope, say along the b axis. We assume again v cx 
= and make use of the above relation between parameters imposed by the experimental anisotropy. The result can 
be considered as an equation giving the hybridization parameter t. Taking T c — 93 K it reads (with t in eV): 



Taking as an example the values we have chosen in the preceding paragraph, together with the BCS isotropic weak 
coupling value A / T c = 1.75, we obtain t = 0.09 eV. However anisotropy as well as strong coupling are expected to 
raise A / T c , and if we rather take the average value found in experiments, namely A / T c « 3 , this raises t up to 
0.16 eV. However we stress again that the low temperature slope is sensitive to the physical parameters in the vicinity 
of the nodes. If we consider that the gap A in the planes has actually some anisotropy, which is quite likely to happen, 
the gap value at the nodes can be lower than the maximum value that we have taken. Similarly we should take into 
account the local value of the mass renormalizations due to interactions, which would be quite hard to obtain. Hence 
it should be clear that we can not escape some appreciable uncertainty in our evaluations. 

Nevertheless the typical value required for the hybridization parameter in order to account for the experiments, 
namely t ss 0.1 eV, is quite satisfactory, taking into account that we have basically no adjustable parameters. Indeed 
this parameter is small enough ( compared to the full band width of the plane or chain bands which is of order 1.6 
eV ) to make our hybridization approach reasonable and consistent. In particular it is small enough for our model to 
explain [jl2| the weak dependence of the critical temperature on impurity content, which is observed experimentally 
and which is one of the major problems of the spin fluctuation model. On the other hand one might worry that the 
hybridization we find is too large, although it is quite obvious from the start that we require a sizeable hybridization, 
since the linear dependence of the penetration depth at low temperature is a fairly strong effect. We will discuss later 
on this point together with our calculations of the penetration depth in the c direction and see that this not the case. 



A y (T)-A y (0) = A y (0) (Av c , y + |A'|v p , y ) 2 
A x (T)-A x (0) A3(0)(Av c , x + |A'|v p , x ) 2 



(13) 




(14) 
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IV. CALCULATIONS OF THE PENETRATION DEPTH OVER THE WHOLE TEMPERATURE RANGE 



In the preceding paragraph we have discussed the low temperature behaviour of the penetration depth and we 
have found that our model can account for experiments in YBCO. The next interesting point is to know if it can 
do the same over the whole temperature range. A priori the answer is likely to be affirmative because the different 
components of this temperature dependence are controlled by different parameters of our model. Thus in principle 
they can be adjusted independently to reproduce the experimental results. Indeed the zero temperature value of the 
penetration depth is an average over the whole Fermi surface of band structure parameters (Fermi velocity, density of 
states), with proper renormalization due to interactions. Therefore it is a global quantity, which in particular does not 
depend on the values of the gaps. On the other hand, as we have already seen explicitely above, the low temperature 
slope of the superfluid density depends only on the band structure parameters and the gaps in the vicinity of the 
nodes, which makes it a local quantity. Finally the slope of p s near T c is a global quantity which depends both on 
band structure parameters and on the gaps on the whole Fermi surface. Since we consider the possibility of strong 
coupling effects, the temperature dependence of the gaps near T c is again an independent variable. The large number 
of parameters is an unpleasant feature of our model, but it is likely to be also an inescapable ingredient of the physics 
of YBa 2 Cu 3 7 . 

Having realized that our model is quite flexible, we can still wonder if we can account for experiment with reasonable 
values of the parameters. In order to explore this problem conveniently we reduce flexibility by making some reasonable 
simplifying assumptions. We take the gaps constant within the plane band and the chain band. The ratio R = |A'|/A 
between the gap A in the plane band and the gap A' in the chain band is one parameter of our calculation. We will 
assume that the temperature dependences of these two gaps are the same and take for it a simple rcnormalized BCS 
form : A(T) = ( A(0) / 1.76 k B T c ) A B cs(T). The ratio A(0) / 1.76 k B T c is another parameter of our calculation. 
For the band structure we take a simple tight binding ansatz which reproduces reasonably well the results obtained 



by band structure calculations pO 21\, and in particular the Fermi surface. We take for the plane band the simple 
tight-binding dispersion relation e& = — 2 t\ ( cos k x + cos k y ) + 2 t% cos k x cos k y — [i , while we use naturally 
e' k = — 2 t[ cos k y — // for the chain band. Actually we find our best overall agreement for t 2 ~ 0. We note that 
the shape of the ( uncoupled ) odd plane band Fermi surface is indeed fairly well described by t 2 ~ 0, except for the 
vicinity of (n, 0) and (0, 7r) which are missed in our simple description ( but since the velocities are small in these 
saddle point regions, this is likely to be unimportant for our purpose ). On the other hand for the even plane band a 
proper description of the Fermi surface requires a larger t 2 . 

Finally, in order to make realistic calculations, we have also included a contribution to the penetration depth 
from this even parity plane band, which is uncoupled to the chains. Wc note that including this contribution makes 
agreement with experiment more difficult to reach. Indeed, since we have no nodes on this piece of the Fermi 
surface, the contribution of this part has an s-wave like temperature dependence and is essentially constant at low 
temperature. This makes it harder to account for the strong slope found experimentally. To perform the calculation, 
we need all the parameters characterizing this band as additional input. These parameters are not related in an 
obvious way to those of the odd parity plane band. Nevertheless we have made a simplifying assumption in order 
to avoid an increase in the number of parameters. The contribution of this plane band to \~ 2 seems isotropic (from 
band structure calculations , its piece of the Fermi surface has essentially a tetragonal symmetry ) . In the a 

direction there is no contribution from the chains in our model, and only the even and odd plane bands contribute. 
We have assumed for simplicity these contributions equal at zero temperature. This means specifically that we have 
taken A~ 2 (0)| el)e „ p i ane = A~ 2 (0)| et)e „ p i ane = A~ 2 (0)| O(W pUne+chain ■ This fixes the T = contribution of the even 
plane band. With respect to its temperature dependence, we take the gap of the even plane equal to the one of the 
odd plane band, namely A(T). These assumptions are quite reasonable, but it is clear that the even plane band is an 
additional source of flexibility of our model, if required. 

Given the large number of parameters still left in our model, we have not proceeded to a systematic exploration in 
the parameter space. We have just made use of our results, obtained in the preceding paragraph, in order to locate 
the region in parameter space for a proper description of the low temperature behaviour, and we have then adjusted 
empirically the parameters to obtain a good averall agreement with experiment. It is possible that better sets of 
parameters can be found. However our purpose in this paper was to check if our model could give a satisfactory 
account for experiment. Since our answer is positive we have not attempted an optimization of our parameters, which 
would be rather meaningless in view of the flexibility mentionned above. Later on we comment on the sensitivity of 
our results to the variation of our parameters. We find a good averall agreement, as shown in Fig. I, with the set of 
parameters h = 0.5 eV, t 2 = - 0.0 eV, t[ = 0.33 eV, n = - 0.58 eV, // = - 0.46 eV , t c = 0.2 eV, |A'|/ A = 0.4 
and 2 A / T c = 2.2 . 
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Let us first consider the T = values of the penetration depths. We find A a = 800 A and \b — 550 A. A priori 
our model is not designed to calculate the absolute values of the penetration depths because it uses a rather primitive 
band structure. It is nevertheless interesting to compare our results to full band structure calculations |29] of the 
plasma frequencies f2 p i which can be related to our penetration depths by Xi — c / Q p i . Their results give A a = 730 
A, A fc = 450 A and A c = 2100 A. We see that our values are rather near these ab initio calculations, which proves that 
our model is quite reasonable. In particular the anisotropy ratio they find is Aj, / A a — 1.62 while we obtain A& / A a 
= 1.45. Both are very near the experimental result 1.6 . This ratio does not vary much when we vary our parameters 
( in a reasonable range). This is easy to understand. This ratio comes merely from the fact that a proper account 
has been made of the chain contribution to the penetration depth along the b axis. This agreement with experiment 
is a strong indication that the chains have to be included in order to get a proper description of the superconducting 
properties of YBCO as we have discussed in the introduction. Nevertheless both our calculation and Ref. |2^] give 
absolute values for Ai which are much lower than the experimental results. The obvious explanation for this is mass 
renormalization, due to the interactions, which we have not taken into account. For an isotropic superconductor 
this would multiply the penetration depth by ( 1 + A ) 1 / 2 , where A is the coupling constant. In the present case 
this renormalization factor is most likely rather anisotropic, but since we have unfortunately no precise idea of this 
anisotropy, let us take it isotropic for the sake of this discussion . If we require for example that it pushes A a to the 
experimental value A a = 1600A, we find = 1100 A, in quite good agreement with experiment. The corresponding 
value of A is 3 which is reasonable for a strong coupling constant. 

Let us turn now to our results for \~ 2 (T) / A~ 2 (0) plotted in Fig.l . As it can be seen we find a very reasonable 
agreement with the low temperature behaviour for both the a and the b directions. However this requires the rather 
high value t c = 0.2 eV which we will discuss later on. The low temperature slopes are naturally sensitive to the other 
band structure parameters because the Fermi velocities in the anticrossing region must be high enough to account 
for experiment. This implies a proper location of the anticrossing region. However, once the proper range for these 
parameters is found, the results are not particularly sensitive to their choice in this range. They are however sensitive 
to the choice of |A'|/ A because this parameter controls the position of the nodes, and the slopes are directly related 
to the orientation and the strength of the Fermi velocities at the nodes. We remark the value |A'|/ A = 0.5 that 
we use is quite reasonable physically : since we expect superconductivity to arise mostly because of interactions in 
the planes, we should find the order parameter smaller in the chains than in the planes. We note also that the 
experimental results displays a linear behaviour over a remarkably large range of temperature since it goes almost 
up to T / T c — 0.5 . A natural explanation for this feature is strong coupling effects which, by increasing 2 A / 
T c , increase the scale for the low temperature behaviour (which is A ) compared to the overall temperature scale 
(which is T c ). This effect is not fully accounted for by our calculations because we have a rather small 2 A / T c — 
2.2 . However a larger value would lower the low temperature slopes and give a disagreement with experiment. This 
situation is due to our simplifying choice of an isotropic gap. Choosing an anisotropic gap would make it possible 
for us to have a large slope ( related to the gap values near the nodes ) and a large linear region (related to the gap 
values farther from the nodes). Finally the agreement is not quite good when we get near T c , but this was expected. 



Indeed experimental results 1 30 1 give a ( 1 - T / T c ) v dependence for A(T) near T c with y w 1/3 , while our mean 
field approch leads us naturally to the standard ( 1 - T / Tc) 1 / 2 dependence. If this experimental result is confirmed 
and the effect is due to a large domain for critical fluctuations ]30| ], including this phyics in our model will naturally 
bring us in better agreement with experiment. 



V. C-AXIS PENETRATION DEPTH 



We come now to the calculation of the c-axis penetration depth. We note first that our simple model is intended to 
describe plane-chain coupling and is at the start a two-dimensional model. In order to make more realistic calculations, 
we have included a k z dependent coupling in its simplest form. However this description is obviously very crude and 
we may fear that it gives a poor account of transport in the c direction. For example our plane-chain coupling does 
not have any (k x ,k y ) dependence. Physically we forget that plane-chain hopping is likely to happen through the 
orbitals of the 04 apical oxygen. Similarly we do not take due account of plane-plane hopping which is obviously 
necessary for transport in the c direction. Clearly a mismatch between plane-plane hopping and plane-chain hopping 
would strongly affect these transport properties ( it has for example been suggested by Xiang and Wheatley |3l| that 
the hopping along the c axis vanishes at the nodes ). It is nevertheless interesting to see what is the result of our 
model for the c-axis penetration depth, and we will see that it does not fare so badly. 

Our raw value at T = is A c = 2250 A. This result compares quite well with the result A c = 2100 A of band structure 
calculations given above. This is quite important because it validates the value of our plane-chain hopping parameter 
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t c = 0.2 eV, which we might have feared to be too large. Therefore this value of t c is in reasonable agreement with 
what comes out of these band structure calculations. Next we have to renormalize this result as we have done above 
for the a and b directions. Taking the same renormalization factor, we obtain A c = 4500 A. This is clearly smaller 
than the experimental result A c 10 000 A although this is not way off the mark. We note that this same problem 
is shared by band structure calculations. 

It is perhaps not too difficult to understand this discrepancy. First the experimental side is probably less secure 
than for the in-plane measurements since A c is not obtained directly, but requires a subtraction of the a or b axis 
contributions which might be difficult to make precisely. Moreover, although experiments agree [l^j32]| to find A c w 
10 000A, they strongly disagree with respect to the temperature dependence. While Mao et al. ]16fl obtain at low 
temperature a linear decrease of A" 2 much stronger than in the a or b direction, Bonn et al. find a much weaker 
dependence than for A a or Af, . On the theoretical side our calculation is performed in the clean limit. While this 
sounds quite reasonable for the in-plane directions, one may reasonably think that transport along the c axis is much 
more sensitive to any kind of scattering or disorder. The importance of incoherent hopping in the c direction has 



been stressed for example by Graf et al. 1 33 1 and Radtke et al. [|34|. Clearly the weakness of coherent hopping in 
the c direction makes incoherent hopping much more relevant than for in-plane transport and this will lead to an 
increase in the penetration depth. Similarly since c axis transport implies electrons going through the chains, it will 
be much more sensitive to chain disorder than in-plane transport. This is confirmed by the strong dependence |35} | 
of the c-axis resistivity on oxygen content in Y BaiCuzO'j-y while the in-plane resisitivity depends much less on y. 
Naturally non stoechiometry brings easily a strong oxygen disorder in the chains, and because of the one dimensional 
nature of these chains, the electronic properties are much more perturbed by this disorder. A related explanation is 
the possible existence of layers of impurities or defects parallel to the planes, which would contribute to block c-axis 
transport. Such a situation would be similar to the one found in films for the ab penetration depth where extrinsic 
"weak links" were shown to increase the penetration depth from the intrinsic value of 1700 A up to 3500 A. The ratio 
between these two values is similar to the discrepancy between our result and the experimental one for A c (0). At the 
same time, these weak links make the low temperature dependence weaker than the intrinsic one. It is tempting to 
attribute the difference between the low temperature behaviour found by Mao et al. and by Bonn et al. to such an 
effect, although this seems in contradiction with the fact that they find essentially the same value for A c (0). 

Let us now consider the temperature dependence that we find for A~ 2 (T)/A~ 2 (0). It can be seen in Fig. 2 that it is 
quite strong. The physical explanation is quite clear. In our model all the k z dependence of the band structure comes 
from hybridization, and hence it is stronger in the anticrossing region. This implies that the largest v^^ z comes from 
this region. Since the nodes of the gap are in this same region, their influence are particularly strong and in particular 
the low temperature slope is large. We notice that our result is very analogous to the one found experimentally |ll| 
by Mao et al. This is naturally satisfactory, but from the above discussion we should be careful not to draw any 
conclusion from this agreement. 

We consider finally a slightly different version of our model which could correct for the discrepancy between our 
result and experiment for A c (0). We may indeed wonder if the large value of A c (0) is not due to a weak plane-plane 
coupling. This possibility has been suggested by Basov et al ||. In order to explore our model in this direction, we 
assume that the plane-plane coupling t p is small. To zeroth order one combination of the plane bands is uncoupled 
and has energy while the orthogonal combination of the plane bands and the chain band are hybridized in the 
usual way with their dispersion relation given by Eq.(8) with tk = \/2 t c . The in- plane penetration depths can be 
calculated to zeroth order. On the other hand the velocity along the c axis is zero at this order, and we have to go 
to first order in t p for the calculation of A c . This gives hvk, z /c = t p sm(k z c) for the even plane band, and hvk, z /c 
= 2 t p t 2 c sin (k z c)/ [ ( e± - e ) 2 + 2 i 2 ] . Therefore we can obtain the experimental value of A c (0) by ajdusting 
t p . We show in Fig. 3 the result of our calculations for the set of parameters t\ = 0.33 eV, t 2 — — 0.0 eV, t[ = 0.4 
eV, (i = - 0.7 eV, // = - 0.56 eV , t c = 0.2 eV, |A'|/ A = 0.5 and 2 A / T c = 3 . As it can be seen the agreement 
with experiment is quite good. For the a direction it is slightly better than in Fig. 1 because the average hybridization 
parameter is larger ( in the present case it is constant while in the preceding case it is proportional to sin (k z c/2) 
) and the ratio A / T c is larger. We note that we find for A C (T) a temperature dependence which is now in good 
agreement with Bonn et al. data p^] . Turning now to the absolute values, our raw results at T — are A a = 1300 A 
and Ab = 650 A. If for example we renormalize Afc to 1000 A (and A a to 2000 A ), we find A c = 10 000 A if we take 
t p — 0.043 eV. This small value justifies our first order expansion. 







VI. CONCLUSION 



In this paper we have investigated the consequences of our coupled plain-chain model of Y BaiCu^Oi for the 
penetration depth and compared them with experiments. We have found that our model accounts quite nicely 
for the anisotropy of the penetration depth in the ab plane at T = (and also for the absolute values, when mass 
renormalization is taken into account). This is obvious qualitatively since the chains contribute only in the b direction. 
But the agreement with experiment is also quite good quantitatively. We reproduce fairly well the linear dependence 
at low temperature and more generally the whole temperature dependence for both the a and the b directions (except 
near T c where experiments do not give a mean field dependence). All this is obtained with a set of parameters 
which arc all quite reasonable physically. In particular we require naturally a sizeable hopping term between plane and 
chain to account for experiment, but the value we use is in agreement with results from band structure calculations. 
Regarding the penetration depth along the c axis, the situation is less satisfactory since two experiments strongly 
disagree with respect to the temperature dependence. We can actually reproduce both temperature dependence by 
making slightly different choices for our model to describe hopping in the c direction. The absolute value of the 
penetration depth disagrees typically by a factor 2 with experiments, but a natural explanation is the important effect 
of impurities and disorder on the transport properties along the c axis. 
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FIG. 1. a) A~ 2 (T) / A" 2 (0) as a function of T/T c . Full line : theory from Eq.(9) with parameters given in the text. Filled 
circles : experimental results from Ref.[17]. b) Same as a) for A^ 2 (T) / A^ 2 (0) as a function of T/T c . 



FIG. 2. A c 2 (T) / A c 2 (0) as a function of T/T c . Full line : theory from Eq.(9) with the same parameters as in Fig.l. Filled 
circles : experimental results from Ref.[32]. 



FIG. 3. a) \a 2 (T) I Aa 2 (0) as a function of T/T c . Full line : theory from Eq.(9) for small plane-plane hopping (see text), 
with parameters given in the text. Filled circles : experimental results from Ref.[17]. b) Same as a) for A^~ 2 (T) / A^ 2 (0) as a 
function of T/T c . c) Same as a) for A~ 2 (T) / A~ 2 (0) as a function of T/T c . Filled circles : experimental results from Ref. [32]. 
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